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TECHNICAL FIELD 

The disclosure below relates to the recovery of face shape from images and 
the facial animation. More particularly, the following description relates to 
determining the head motion between two images based on symmetrical features 
in the two images. 

BACKGROUND 

One of the most interesting and difficult problems in computer graphics is 
the effortless generation of realistic looking, animated human face models. 
Animated face models are essential to computer games, film making, online chat, 
virtual presence, video conferencing, etc. So far, the most popular commercially 
available tools have utilized laser scanners. Not only are these scanners expensive, 
the data are usually quite noisy, requiring hand touchup and manual registration 
prior to animating the model. Because inexpensive computers and cameras are 
widely available, there is a great interest in producing face models directly from 
images. In spite of progress toward this goal, the available techniques are either 
manually intensive or computationally expensive. 

Facial modeling and animation has been a computer graphics research topic 
for over 25 years [6, 16, 17, 18, 19, 20, 21, 22, 23, 27, 30, 31, 33]. The reader is 
referred to Parke and Waters' book [23] for a complete overview. 

Lee et al. [17, 18] developed techniques to clean up and register data 
generated from laser scanners. The obtained model is then animated using a 
physically based approach. 

DeCarlo et al. [5] proposed a method to generate face models based on face 
measurements randomly generated according to anthropometric statistics. They 
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showed that they were able to generate a variety of face geometries using these 
face measurements as constraints. 

A number of researchers have proposed to create face models from two 
views [1, 13, 4]. They all require two cameras which must be carefully set up so 
that their directions are orthogonal. Zheng [37] developed a system to construct 
geometrical object models from image contours, but it requires a turn-table setup. 

Pighin et al. [26] developed a system to allow a user to manually specify 
correspondences across multiple images, and use vision techniques to computer 3D 
reconstructions. A 3D mesh model is then fit to the reconstructed 3D points. They 
were able to generate highly realistic face models, but with a manually intensive 
procedure. 

Blanz and Vetter [3] demonstrated that linear classes of face geometries and 
images are very powerful in generating convincing 3D human face models from 
images. Blanz and Vetter used a large image database to cover every skin type. 

Kang et al. [14] also use linear spaces of geometrical models to construct 
3D face models from multiple images. But their approach requires manually 
aligning the generic mesh to one of the images, which is in general a tedious task 
for an average user. 

Fua et al. [8] deform a generic face model to fit dense stereo data, but their 
face model contains a lot more parameters to estimate because basically all of the 
vertexes are independent parameters, plus reliable dense stereo data are in general 
difficult to obtain with a single camera. Their method usually takes 30 minutes to 
an hour, while ours takes 2-3 minutes. 

Guenter et al. [9] developed a facial animation capturing system to capture 
both the 3D geometry and texture image of each frame and reproduce high quality 
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facial animations. The problem they solved is different from what is addressed 
here in that they assumed the person's 3D model was available and the goal was to 
track the subsequent facial deformations. 

SUMMARY 

The system described below provides improved procedures to estimate head 
motion between two images of a face. A procedure is described that first, 
identifies locations of a number of distinct facial features in two images. The 
procedure of estimating head motion with respect to these locations corresponds to 
the determination of a number of unknown quantities. For example, these 
identified locations can correspond to the eye corners, mouth corners and nose tip. 

Next, the locations are converted into as a set of physical face parameters 
based on the symmetry of the identified distinct facial features. The set of physical 
parameters reduces the number of unknowns as compared to the number of 
equations used to determine the unknowns. This reduction in unknowns and 
relative increase in the number of equations used to determine the unknowns 
increases the redundancy and thus, the robustness of the head motion estimation. 

Finally, the points corresponding to the identified points between the two 
images are used to solve the head motion. The head motion estimation operation 
includes: (a) estimating each of the set of physical parameters, (b) estimating a first 
head pose transform corresponding to the first image, and (c) estimating a second 
head pose transform corresponding to the second image. One of the physical 
parameters is set to a constant due to the fact that a scale cannot be determined 
from the two images. 
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Optionally, an inequality constraint can be placed on a particular physical 
parameter, such that the parameter is constrained within a predetermined minimum 
and maximum value. For example, one description provides that the parameter 
corresponds to a nose tip. The predetermined minimum value is zero (0) and the 
predetermined maximum value is a reasonable value based on absolute values of 
other locations - recall that each location identifies a facial feature. The inequality 
constraint is converted to an equality constraint by using a penalty function. Then, 
the inequality constraint is used during the head motion estimation to add 
additional robustness to the motion estimation. 

In yet another description, a procedure is provided to use the head motion 
estimation described above as an initial estimation in combination with a feature 
matching algorithm. To accomplish this, the procedure first involves identifying 
locations of a plurality of distinct facial features in the two images, the locations 
corresponding to a number of unknowns determined upon estimation of head 
motion. For example, these identified locations can correspond to the eye corners, 
mouth corners and nose tip. 

Next, the identified locations are converted into a set of physical face 
parameters based on the symmetry of the identified distinct facial features, the set 
of physical parameters reducing the number of unknowns. Next, a first set of 
matched points (corresponding to previously identified distinct facial features) is 
used to determine the head motion and the physical facial parameters. Finally, a 
second set of matched points from a feature matching algorithm is incorporated, 
together with the first set of matched points, to refine the estimation of the head 
motion and the physical facial parameters. These novel procedures provide 
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additional robustness to systems that estimate head motion during animation 
modeling. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a block diagram of a computer system capable of performing the 
operations described below. 

Fig. 2 illustrates how to mark facial features on an image. 

Figs. 3, 5, 6 are flow charts showing sequences of actions for creating a 3D 
face model. 

Fig. 4 shows the selection of different head regions as described below. 
Fig. 7 illustrates a coordinate system that is based on symmetry between 
selected feature points on an image. 

DETAILED DESCRIPTION 

The following description sets forth a specific embodiment of a 3D 
modeling system that incorporates elements recited in the appended claims. The 
embodiment is described with specificity in order to meet statutory requirements. 
However, the description itself is not intended to limit the scope of this patent. 
Rather, the claimed invention might eventually be embodied in other ways, to 
include different elements or combinations of elements similar to the ones 
described in this document, in conjunction with other present or future 
technologies. 
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System Overview 

Fig. 1 shows components of our system. The equipment includes a 
computer 10 and a video camera 12. The computer is a typical desktop, laptop, or 
similar computer having various typical components such as a keyboard/mouse, 
display, processor, peripherals, and computer-readable media on which an 
operating system and application programs are stored and from which the 
operating system and application programs are executed. Such computer-readable 
media might include removable storage media, such as floppy disks, CDROMs, 
tape storage media, etc. The application programs in this example include a 
graphics program designed to perform the various techniques and actions 
described below. 

The video camera is an inexpensive model such as many that are widely 
available for Internet videoconferencing. We assume the intrinsic camera 
parameters have been calibrated, a reasonable assumption given the simplicity of 
calibration procedures [36]. 

Data Capture 

The first stage is data capture. The user takes two images with a small 
relative head motion, and two video sequences: one with the head turning to each 
side. Alternatively, the user can simply turn his/her head from left all the way to 
the right, or vice versa. In that case, the user needs to select one approximately 
frontal view while the system automatically selects the second image and divides 
the video into two sequences. In the seque, we call the two images the base 
images. 
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The user then locates five markers in each of the two base images. As 
shown in Fig. 2, the five markers correspond to the two inner eye corners 20, nose 
tip 21, and two mouth corners 22. 

The next processing stage computes the face mesh geometry and the head 
pose with respect to the camera frame using the two base images and markers as 
input. 

The final stage determines the head motions in the video sequences, and 
blends the images to generate a facial texture map. 

Notation 

We denote the homogeneous coordinates of a vector x by x, i.e., the 
homogeneous coordinates of an image point m = (u,v) T are m = (u,v,l f , and those 
of a 3D point p =(x,y,zf are p = (x,y,z,l) T . A camera is described by a pinhole 
model, and a 3D point p and its image point m are related by 

Am = APQp 

where A is a scale, and A, P, and Q are given by 
A = 
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The elements of matrix A are the intrinsic parameters of the camera and matrix A 
maps the normalized image coordinates to the pixel image coordinates (see e.g. 
[7]). Matrix P is the perspective projection matrix. Matrix q is the 3D rigid 
transformation (rotation R and translation t) from the object/world coordinate 
system to the camera coordinate system. When two images are concerned, a prime 
' is added to denote the quantities related to the second image. 
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The fundamental geometric constraint between two images is known as the 
epipolar constraint [7, 35]. It states that in order for a point m in one image and a 
point m' in the other image to be the projections of a single physical point in space, 
or in other words, in order for them to be matched, they must satisfy 

m' r A' r EA- 1 rn = 0 

where E = [t r ] x R r is known as the essential matrix, (R r tJ is the relative motion 
between the two images, and [t,] x is a skew symmetric matrix such that 
t r x v = [t r ] x v for any 3D vector v. 

Linear Class of Face Geometries 

Instead of representing a face as a linear combination of real faces or face 
models, we represent it as a linear combination of a neutral face model and some 
number of face metrics, where a metric is a deformation vector that linearly 
deforms a face in a certain way, such as to make the head wider, make the nose 
bigger, etc. Each deformation vector specifies a plurality of displacements 
corresponding respectively to the plurality of 3D points of the neutral face model. 

To be more precise, let's denote the face geometry by a vector 
s = ( v i r >-<) r , where v, =(X i ,Y i ,Z l ) T , (z = l,...,«) are the vertices, and a metric by 
a vector M = {S\ x ,...,Sv n ) T , where Sv i =(SX i ,SY i ,SZ l ) T . Given a neutral face 
S° =(\° r ,...y 0 n T ) r , and a set of m metrics M J =(Sv{ T , the linear space of 
face geometries spanned by these metrics is 

m 

S = S° +^ Cj M J subject to c y e [/,.,«,] 

where c/s are the metric coefficients and lj and uj are the valid range of c } . 
In our implementation, the neutral face and all the metrics are designed by an 
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artist, and it is done only once. The neutral face contains 194 vertices and 360 
triangles. There are 65 metrics. 

Image Matching and 3D Reconstruction 

We now describe our techniques to determine the face geometry from just 
two views. The two base images are taken in a normal room by a static camera 
while the head is moving in front. There is no control on the head motion, and the 
motion is unknown. We have to determine first the motion of the head and match 
some pixels across the two views before we can fit an animated face model to the 
images. However, some preprocessing of the images is necessary. 

Determining Facial Portions of the Images 

Fig. 3 shows actions performed to distinguish a face in the two selected 
images from other portions of the images. 

There are at least three major groups of objects undergoing different 
motions between the two views: background, head, and other parts of the body 
such as the shoulder. If we do not separate them, there is no way to determine a 
meaningful head motion, since the camera is static, we can expect to remove the 
background by subtracting one image from the other. However, as the face color 
changes smoothly, a portion of the face may be marked as background. Another 
problem with the image subtraction technique is that the moving body and the head 
cannot be distinguished. 

An initial step 100 comprises using image subtraction to create the first 
mask image in which pixels having different colors in the two base images are 
marked. 
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A step 101 comprises identifying locations of a plurality of distinct facial 
features in the base images. In this example, the user does this manually, by 
marking the eyes, nose, and mouth, as described above and shown in Fig. 2. 
Automated techniques could also be used to identify these points. 

A step 102 comprises calculating a range of skin colors by sampling the 
base images at the predicted portions, or at locations that are specified relative to 
the user-indicated locations of the facial features. This allows us to build a color 
model of the face skin. We select pixels below the eyes and above the mouth, and 
computer a Gaussian distribution of their colors in the RGB space. If the color of 
a pixel matches this face skin color model, the pixel is marked as a part of the face. 

A step 103 comprises creating a second mask image that marks any image 
pixels having colors corresponding to the calculated one or more skin colors. 

Either union or intersection of the two mask images is not enough to locate 
the face because it will include either too many (e.g., including undesired moving 
body) or too few (e.g., missing desired eyes and mouth) pixels. Since we already 
have information about the position of eye corners and mouth corners, we initially 
predict the approximate boundaries of the facial portion of each image, based on 
the locations identified by the user. More specifically, step 104 comprises 
predicting an inner area and an outer area of the image. The outer area 
corresponds roughly to the position of the subject's head in the image, while the 
inner area corresponds roughly to the facial portion of the head. 

Fig. 4 shows these areas, which are defined as ellipses. The inner ellipse 23 
covers most of the face, while the outer ellipse 24 is usually large enough to 
enclose the whole head. Let d e be the image distance between the two inner eye 
corners, and d em , the vertical distance between the eyes and the mouth. The width 
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and height of the inner ellipse are set to 5d e and 3d em . The outer ellipse is 25% 
larger than the inner one. 

In addition, step 104 includes predicting or defining a lower area of the 
image that corresponds to a chin portion of the head. The lower area aims at 
removing the moving body, and is defined to be 0.6d em below the mouth. 

Within the inner ellipse, a "union" or "joining" operation 105 is used: we 
note all marked pixels in the first mask image and also any unmarked pixels of the 
first mask image that correspond in location to marked pixels in the second mask 
image. Between the inner and outer ellipses (except for the lower region),an 
image subtraction operation 106 is used: we note which pixels (marked or 
unmarked) of the image have different colors relative to correspondingly located 
pixels in the other image. In the lower part, we use an "intersection" operation 
107: we note any marked pixels in the first mask image that correspond in location 
to marked pixels in the second mask image. 

A step 108 comprises forming a final mask image that marks the noted 
pixels as being part of the head. This involves joining the mask image to the mask 
image of the other image. More specifically, the corresponding noted pixels of 
each base image are logically OR's to create the final maskimage. 

Corner Matching and Motion Determination 

One popular technique of image registration is optical flow [12, 2], which is 
based on the assumption that the intensity/color is conserved. This is not the case 
in our situation: the color of the same physical point appears to be different in 
images because the illumination changes when the head is moving. We therefore 
resort to a feature-based approach that is more robust to intensity/color variations. 
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It consists of the following steps: (i) detecting comers in each image; (ii) matching 
corners between the two images; (iii) detecting false matches based on a robust 
estimation technique; (iv) determining the head motion; (v) reconstructing matched 
points in 3D space. 

Fig. 5 shows the sequence of operations. 

Corner Detection. In a step 110, we use the Plessey corner detector, a 
well-known technique in computer vision [10]. It locates corners corresponding to 
high curvature points in the intensity surface if we view an image as a 3D surface 
with the third dimension being the intensity. Only corners whose pixels are white 
in the mask image are considered. 

Corner Matching. In a step 111, for each corner in the first image we 
choose an 11x11 window centered on it, and compare the window with windows 
of the same size, centered on the corners in the second image. A zero-mean 
normalized cross correlation between two windows is computed [7]. If we 
rearrange the pixels in each window as a vector, the correlation score is equivalent 
to the cosine angle between two intensity vectors. It ranges from -1, for two 
windows which are not similar at all, to 1, for two windows which are identical. If 
the largest correlation score exceeds a prefixed threshold (0.866 in our case), then 
that corner in the second image is considered to be the match candidate of the 
corner in the first image. The match candidate is retained as a match if and only if 
its match candidate in the first image happens to be the corner being considered. 
This symmetric test reduces many potential matching errors. 

False Match Detection. Operation 112 comprises detecting and discarding 
false matches. The set of matches established so far usually contains false matches 
because correlation is only a heuristic. The only geometric constraint between two 
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images is the epipolar constraint m' r A'" r EA _1 m = 0 . If two points are correctly 
matched, they must satisfy this constraint, which is unknown in our case. 
Inaccurate location of corners because of intensity variation of lack of string 
texture features is another source of error. In a step 109, we use the technique 
described in [35] to detect both false matches and poorly located corners and 
simultaneously estimate the epipolar geometry (in terms of the essential matrix E). 
That technique is based on a robust estimation technique known as the least 
median squares [28], which searches in the parameter space to find the parameters 
yielding the smallest value for the median of squared residuals computer for the 
entire data set. Consequently, it is able to detect false matches in as many as 
49.9% of the whole set of matches. 

Motion Estimation 

In a step 113, we compute an initial estimate of the relative head motion 
between two images, denoted by rotation R r and translation t,. If the image 
locations of the identified feature points are precise, one could use a five-point 
algorithm to compute camera motion from Matrix E [7, 34]. Motion (R r , t r ) is then 
re-estimated with a nonlinear least-squares technique using all remaining matches 
after having discarded the false matches [34]. 

However, the image locations of the feature point are not usually precise. A 
human typically cannot mark the feature points with high precision. An automatic 
facial feature detection algorithm may not produce perfect results. When there are 
errors, a five-point algorithm is not robust even when refined with a well-known 
bundle adjustment technique. 
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For each of the five feature points, its 3D coordinates (x, y, z) coordinates 
need to be determined - fifteen (1 5) unknowns. Then, motion vector (R r , t,) needs 
to be determined - adding six (6) more unknowns. One unknown quantity is the 
magnitude, or global scale, which will never be determined from images alone. 
Thus, the number of unknown quantities that needs to be determined is twenty 
(i.e., 15+6-1 = 20). The calculation of so many unknowns further reduces the 
robustness of the five point-tracking algorithm. 

To substantially increase the robustness of the five point algorithm, a new 
set of parameters is created. These parameters take into consideration physical 
properties of the feature points. The property of symmetry is used to reduce the 
number of unknowns. Additionally, reasonable lower and upper bounds are placed 
on nose height and are represented as inequality constraints. As a result, the 
algorithm becomes more robust. Using these techniques, the number of unknowns 
is significantly reduced below 20. 

Even though the following algorithm is described with respect to five 
feature points, it is straightforward to extend the idea to any number of feature 
points less than or greater than five feature points for improved robustness. 
Additionally, the algorithm can be applied to other objects besides a face as long as 
the other objects represent some level of symmetry. Head motion estimation is 
first described with respect to five feature points. Next, the algorithm is extended 
to incorporate other image point matches obtained from image registration 
methods. 

Head Motion Estimation from Five Feature Points. Fig. 7 illustrates the 
new coordinate system used to represent feature points. E { 202, E 2 204, M r 206, 
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M 2 208, and N 210 denote the left eye comer, right eye corner, left mouth corner, 
right mouth corner, and nose top, respectively. A new point E 212 denotes the 
midpoint between eye corners E l5 E 2 and a new point M 214 identifies the 
midpoint between mouth corners M h M 2 . Notice that human faces exhibit some 
strong structural properties. For example, the left and right sides of a human face 
are very close to being symmetrical about the nose. Eye corners and mouth 
corners are almost coplanar. Based on these symmetrical characteristics, the 
following reasonable assumptions are made: 

(1) A line E^ connecting the eye corners Ei and E 2 is parallel to a line 
M]M 2 connecting the mouth corners. 

(2) A line centered on the nose (e.g., line EOM when viewed straight on or 
lines NM or NE when viewed from an angle as shown) is perpendicular 
to mouth line MxM 2 and to eye line E^. 

Let n be the plane defined by E l5 E2, Mi and M 2 . Let O 216 denote the 
projection of point N on plane n . Let Q 0 denote the coordinate system, which is 
originated at 0 with ON as the z-axis, OE as the y -axis; the x-axis is defined 
according to the right-hand system. In this coordinate system, based on the 
assumptions mentioned earlier, we can define the coordinates of Ei, E 2 , M u M 2 , N 
as (-a, b, 0) T , (a, b, 0) T , (-d, -c, 0) T , (d, -c, 0) T , (0, 0, e) T , respectively. 

By redefining the coordinate system, the number of parameters used to 
define five feature points is reduced from nine (9) parameters for generic five 
points to five (5) parameters for five feature points in this local coordinate system. 

Let t denote the coordinates of O under the camera coordinate system, and 
R the rotation matrix whose three columns are vectors of the three coordinate axis 
of Oq. For each point p e {E u E 2 ,M b M 2 , N}, its coordinate under the camera 
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coordinate system is Rp +t. We call (R, t) the head pose transform. Given two 
images of the head under two different poses (assume the camera is static), let 
(R, t) and (R', t') be their head pose transforms. For each point pj g { E h E 2 ,M l5 
M 2 , N}, if we denote its image point in the first view by m; and that in the second 
view by m\, we have the following equations: 
prq/-(Rp i + t) = m i (1) 

and 

proj(R>p i + V) = m\ (2) 
where pro] is the perspective projection. Notice that we can fix one of the 
coordinates a, b, c, d, since the scale of the head size cannot be determined from 
the images. As is well known, each pose has six (6) degrees of freedom. 
Therefore, the total number of unknowns is sixteen (16), and the total number of 
equations is 20. If we instead use their 3D coordinates as unknowns as in any 
typical bundle adjustment algorithms, we would end up with 20 unknowns and 
have the same number of equations. By using the generic properties of the face 
structure, the system becomes over-constrained, making the pose determination 
more robust. 

To make the system even more robust, we add an inequality constraint on e. 
The idea is to force e to be positive and not too large compared to a, b, c, d. In the 
context of the face, the nose is always out of plane n. In particular, we use the 
following inequality: 

0 < e < 3a (3) 
Three (3) is selected as the upper bound of el a simply because it seems reasonable 
and it works well. The inequality constraint is finally converted to equality 
constraint by using a penalty function. 
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e*e ife<0 

0 if 0 < e < 3a ^ 
(e-3a)*(e-3a) ife>3a 



In summary, based on equations (1), (2) and (4), we estimate a, b, c 9 d, e, 
(R, t) and (R' ? t') by minimizing 

5 

=E W '(K -proj<$Q t +t)f +|K-^(R'p i +t'f) + w n P nose (5) 

where wfs and w n are the weighting factors, reflecting the contribution of each 
term. In our case, w t = 1 except for the nose term which has a weight of 0.5 
because it is usually more difficult to locate the nose top than other feature points. 
The weight for penalty w n is set to 10. The objective function (5) is minimized 
using a Levenberg-Marquardt method [40]. More precisely, as mentioned earlier, 
we set a to a constant during minimization since the global head size cannot be 
determined from images. 

Incorporating Image Point Matches. If we estimate camera motion using 
only the five user marked points, the result is sometimes not very accurate because 
the markers contain human errors. In this section, we describe how to incorporate 
the image point matches (obtained by any feature matching algorithm) to improve 
precision. 

Let (nij, m'j) (j = 1 ... K) be the K point matches, each corresponding to the 
projections of a 3D point pj according to the perspective projection (1) and (2). 3D 
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points pj's are unknown, so they are estimated. Assuming that each image point is 
extracted with the same accuracy, we can estimate a, b, c, d, e, (R, t), (R\ t'), and 
{ Pj } 0 = 1 - K ) by minimizing 

K 

F = Fs P ts +^X<J m y -P r °J( R Pj +*)| +Wj-proj(R'pj + t')f) (6) 

where F 5pts is given by (5), and w v is the weighting factor. We set w p = 1 by 
assuming that the extracted points have the same accuracy as those of eye corners 
and mouth corners. The minimization can again be performed using a Levenberg- 
10 Marquardt method. This is a quite large minimization problem since we need to 
n estimate 16 +3 K unknowns, and therefore it is computationally quite expensive 

12 especially for large K. Fortunately as shown in [37], we can eliminate the 3D 

13 points using a first order approximation. The following term 

14 

15 | ||m y ~proj'(Rp j +t)|p +\m , J -proj(R'p. +t')|f 

16 

l? I can be shown to be equal, under the first order approximation, to 



18 



19 



20 



21 



22 



(mj Em ) 2 



mj E r ZZ r Em y +mj E r ZZ r Em; 



where mi = [m y r ,lf,m; = [my ,lf,Z = 
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, and E is the essential matrix to be 
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Let (R r , tr) be the relative motion between two views. It is easy to see that 
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R r = R'R T , and 
tr = V -R'R T t. 

Furthermore, let's define a 3 x3 antisymmetric matrix [tj* such that [tj x x = t, x x 
for any 3D vector x. The essential matrix is then given by 

E = [t r ]xR r (7) 
which describes the epipolar geometry between two views [7]. 
In summary, the objective function (6) becomes 

F = F 5pt5 + -i± _ (8) 

>» mj E r ZZ r Em, + mj E r ZZ r Em; 

Notice that this is a much smaller minimization problem. We only need to 
estimate 16 parameters as in the five-point problem (5), instead of 16 +3 K 
unknowns. 

To obtain a good initial estimate, we first use only the five feature points to 
estimate the head motion by using the algorithm described in Section 2. Thus we 
have the following two step algorithm: 

Stepl. Set w p =0. Solve minimization problem 8. 

Step2. Set w p =l. Use the results of stepl as the initial estimates. Solve 
minimization problem (8). 

Notice that we can apply this idea to the more general cases where the 
number of feature points is not five. For example, if there are only two eye corners 
and mouth corners, we'll end up with 14 unknowns and 16 +3 K equations. Other 
symmetric feature points (such as the outside eye corners, nostrils, and the like) 
can be added into equation 8 in a similar way by using the local coordinate system 
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Head Motion Estimation Results 

In this section, we show some test results to compare the new algorithm 
with the traditional algorithms. Since there are multiple traditional algorithms, we 
chose to implement the algorithm as described in [34]. It works by first computing 
an initial estimate of the head motion from the essential matrix [7], and then re- 
estimate the motion with a nonlinear least-squares technique. 

We have run both the traditional algorithm and the new algorithm on many 
real examples. We found many cases where the traditional algorithm fails while the 
new algorithm successfully results in reasonable camera motions. When the 
traditional algorithm fails, the computed motion is completely bogus, and the 3D 
reconstructions give meaningless results. But the new algorithm gives a 
reasonable result. We generate 3D reconstructions based on the estimated motion, 
and perform Delauney triangulation. 

We have also performed experiments on artificially generated data. We 
arbitrarily select 80 vertices from a 3D face model and project its vertices on two 
views (the head motion is eight degrees apart). The image size is 640 by 480 
pixels. We also project the five 3D feature points (eye corners, nose top, and mouth 
corners) to generate the image coordinates of the markers. We then add random 
noises to the coordinates (u, v) of both the image points and the markers. The 
noises are generated by a pseudo-random generator subject to Gausian distri- 
bution with zero mean and variance ranging from 0.4 to 1.2. We add noise to the 
marker's co-ordinates as well. The results are plotted in Figure 3. The blue curve 
shows the results of the traditional algorithm and the red curve shows the results of 
our new algorithm. The horizontal axis is the variance of the noise distribution. 
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The vertical axis is the difference between the estimated motion and the actual 
motion. The translation vector of the estimated motion is scaled so that its 
magnitude is the same as the actual motion. The difference between two rotations 
is measured as the Euclidean distance between the two rotational matrices. 

We can see that as the noise increases, the error of the traditional algorithm 
has a sudden jump at certain point. But, the errors of our new algorithm grow 
much more slowly. 

3D Reconstruction . 

Fig. 6 illustrates a step 114, where matched points are reconstructed in 3D 
space with respect to the camera frame at the time when the first base image was 
taken. Let (m, m') be a couple of matched points, and p be their corresponding 
point in space. 3D point p is estimated such that |m-m| 2 +||m'-m'|| 2 is 
minimized, where m and hY are projections of p in both images according to the 
equation Am = APQp . 

3D positions of the markers are determined in the same way. 

Fitting a Face Model 

This stage of processing creates a 3D model of the face. The face model 
fitting process consists of two steps: fitting to 3D reconstructed points and fine 
adjustment using image information. 

3D Fitting 

A step 120 comprises constructing a realistic 3D face model from the 
reconstructed 3D image calculated in step 111. Given a set of reconstructed 3D 
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points from matched comers and markers, the fitting process applies a combination 

of deformation vectors to a pre-specified, neutral face model, to deform the neutral 

face model approximately to the reconstructed face model. The technique searches 

for both the pose of the face and the metric coefficients to minimize the distances 

from the reconstructed 3D points to the neutral face mesh. The pose of the face is 

fsR i) 

the transformation T= from the coordinate frame of the neutral face mesh 

v° v 

to the camera frame, where R is a 3x3 rotation matrix, t is a translation, and s is a 
global scale. For any 3D vector p, we use notation T(p) = sRp + 1 . 

The vertex coordinates of the face mesh in the camera frame is a function of 
both the metric coefficients and the pose of the face. Given metric coefficients (c,, 
. ■ .,c m ) and pose T, the face geometry in the camera frame is given by 

f n \ 

S = T 



Since the face mesh is a triangular mesh, any point on a triangle is a linear 
combination of the three triangle vertexes in terms of barycentric coordinates. So 
any point on a triangle is also a function of T and metric coefficients. 
Furthermore, when T is fixed, it is simply a linear function of the metric 
coefficients. 

Let (Pi, P2, • • , Pk) be the reconstructed corner points, and (q b q 2 , . . . , q 5 ) 
be the reconstructed markers. Denote the distance from p ; to the face mesh S by 
d(p i} S). Assume marker q y corresponds to vertex v m . of the face mesh, and 
denote the distance between q, and \ m . by rf(q y , ) . The fitting process consists 
of finding pose T and metric coefficients {c x , . . . , c n } by minimizing 

.=1 y=i J 

where w t is a weighting factor. 
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To solve this problem, we use an iterative closest point approach. At each 
iteration, we first fix T. For each p„ we find the closest point g ; on the current face 
mesh S. We then minimize £ V 2 (/>,•>£) + X^ta,-,^) . We set w t to be 1 at 
the first iteration and 1.0/l+d 2 (p h g,)) in the subsequent iterations. The reason for 
using weights is that the reconstruction from images is noisy and such a weight 
scheme is an effective way to avoid overfitting to the noisy data [8]. Since both g ; 
and v mj are linear functions of the metric coefficients for fixed T, the above 
problem is a linear least square problem. We then fix the metric coefficients, and 
solve for the pose. To do that, we recompute g ; using the new metric coefficients. 
Given a set of 3D corresponding points (p fj & ) and (q . v m ), there are well known 
algorithms to solve for the pose. We use the quaternion-based technique described 
in [1 1]. To initialize this iterative process, we first use the 5 markers to compute 
an initial estimate of the pose. In addition, to get a reasonable estimate of the head 
size, we solve for the head-size related metric coefficients such that the resulting 
face mesh matches the bounding box of the reconstructed 3D points. 
Occasionally, the corner matching algorithm may produce points not on the face. 
In that case, the metric coefficients will be out of the valid ranges, and we throw 
away the point that is the most distant from the center of the face. We repeat this 
process until metric coefficients become valid. 

Fine Adjustment Using Image Information 

After the geometric fitting process, we have now a face mesh that is a close 
approximation to the real face. To further improve the result, we perform a search 
130 for silhouettes and other face features in the images and use them to refine the 
face geometry. The general problem of locating silhouettes and face features in 
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images is difficult, and is still a very active research area in computer vision. 
However, the face mesh that we have obtained provides a good estimate of the 
locations of the face features, so we only need to perform search in a small region. 

We use the snake approach [15] to computer the silhouettes of the face. The 
silhouette of the current face mesh is used as the initial estimate. For each point on 
this piecewise linear curve, we find the maximum gradient location along the 
normal direction within a small range (10 pixels each side in our implementation). 
Then we solve for the vertexes (acting as control points) to minimize the total 
distance between all the points and their corresponding maximum gradient 
locations. 

We use a similar approach to find the upper lips. 

To find the outer eye corner (not marked), we rotate the current estimate of 
that eye corner (given by the face mesh) around the marked eye corner by a small 
angle, and look for the eye boundary using image gradient information. This is 
repeated for several angles, and the boundary point that is the most distant to the 
marked corner is chosen as the outer eye corner. 

We could also use the snake approach to search for eyebrows. However, 
our current implementation uses a slightly different approach. Instead of 
maximizing image gradients across contours, we minimize the average intensity of 
the image area that is covered by the eyebrow triangles. Again, the vertices of the 
eyebrows are only allowed to move in a small region bounded by their neighboring 
vertices. This has worked very robustly in our experiments. 

We then use the face features and the image silhouettes as constraints in our 
system to further improve the mesh, in a step 131. Notice that each vertex on the 
mesh silhouette corresponds to a vertex on the image silhouette. We cast a ray 
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from the camera center through the vertex on the image silhouette. The projection 
of the corresponding mesh vertex on this ray acts as the target position of the mesh 
vertex. Let v be the mesh vertex and h the projection. We have equation v = h. 
For each face feature, we obtain an equation in a similar way. These equations are 
added to equation (5). The total set of equations is solved as before, i.e., we first 
fix the post T and use a linear least square approach to solve the metric 
coefficients, and then fix the metric coefficients while solving for the pose. 

Face Texture From Video Sequence 

Now we have the geometry of the face from only two views that are close to 
the frontal position. For the sides of the face, the texture from the two images is 
therefore quite poor or even not available at all. Since each image only covers a 
portion of the face, we need to combine all the images in the video sequence to 
obtain a complete texture map. This is done by first determining the head pose for 
the images in the video sequence and then blending them to create a complete 
texture map. 

Determining Head Motions in Video Sequences 

Fig. 6 shows operations in creating a texture map. In an operation 140, 
successive images are first matched using the same corner detection, corner 
matching, and false match detection techniques described above. We could 
combine the resulting motions incrementally to determine the head pose. 
However, this estimation is quite noisy because it is computed only from 2D 
points. As we already have the 3D face geometry, a more reliable pose estimation 
can be obtained by combining both 3D and 2D information, as follows. 
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In an operation 141, the pose of each successive image is determined. Let 

us denote the first base image by I 0 . This base image comprises one of the two 

initial still images, for which the pose is already known. Because we know the 

pose of the base image, we can determine the 3D position of each point in the base 

image relative to the facial model that has already been computed. 

We will denote the images on the video sequences by I h . . . , /„. The 

(R t ^ 

relative head motion from I iA to I t is given by R = * " , and the head pose 

v° U 

corresponding to image I t with respect to the camera frame is denoted by Q,-. The 
technique works incrementally, starting with I 0 and I h For each pair of images (I u 
1, I t ), we perform a matching operation to match points of image with 
corresponding points in I iA . This operation uses the corner matching algorithm 
described above. We then perform a minimization operation, which calculates the 
pose of Ii such that projections of 3D positions of the matched points of I iA onto I t 
coincide approximately with the corresponding matched points of I t . More 
specifically, the minimization operation minimizes differences between the 
projections of 3D positions of the matched points of I lA onto I t and the 
corresponding matched points of I t . Let us denote the matched corner pairs as 
t m j' m 'j) | .7 = 1,— ,/}- For each nij in I iA , we cast a ray from the camera center 
through mj, and compute the intersection Xj of that ray with the face mesh 
corresponding to image According to the equation Am = APQp , R ; is subject 
to the following equations 

APR t x.=A..m' J fory=l, ...,/ 
where A, P, Xj and m^. are known. Each of the above equations gives two 
constraints on R,-. We compute R ; with a technique described in [37], which 
minimizes the sum of differences between each pair of matched points [m j ,rri: j ). 
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After R y is computed, the head pose for image 7, in the camera frame is given by Q t 
= Rj The head pose Q 0 is known from previous calculations involving the 
two still images. 

In general, it is inefficient to use all the images in the video sequence for 
texture blending, because head motion between two consecutive frames is usually 
very small. To avoid unnecessary computation, the following process is used to 
automatically select images from the video sequence. Let us call the amount of 
rotation of the head between two consecutive frames the rotation speed. Ifs is the 
current rotation speed and a is the desired angle between each pair of selected 
images, the next image is selected a/s frames away. In our implementation, the 
initial guess of the rotation speed is set to 1 degree/frame and the desired 
separation angle is equal to 5 degrees. 

Texture Blending 

Operation 142 is a texture blending operation. After the head pose of an 
image is computed, we use an approach similar to Pighin et al.'s method [26] to 
generate a view independent texture map. We also construct the texture map on a 
virtual cylinder enclosing the face model. But instead of casting a ray from each 
pixel to the face mesh and computing the texture blending weights on a pixel by 
pixel basis, we use a more efficient approach. For each vertex on the face mesh, 
we computed the blending weight for each image based on the angle between 
surface normal and the camera direction [26]. If the vertex is invisible, its weight 
is set to 0.0. The weights are then normalized so that the sum of the weights over 
all the images is equal to 1.0. We then set the colors of the vertexes to be their 
weights, and use the rendered image of the cylindrical mapped mesh as the weight 
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map. For each image, we also generate a cylindrical texture map by rendering the 
cylindrical mapped mesh with the current image as texture map. Let C ( and 
Wi (I = 1, . . . , k) be the cylindrical texture maps and the weight maps. Let D be 
the final blended texture map. For each pixel (u, v), its color on the final blended 
texture map is 

C(u,v) = f j W i (u,v)C i (u,v). 

i=l 

Because the rendering operations can be done using graphics hardware, this 
approach is very fast. 

User Interface 

We have built a user interface to guide the user through collecting the 
required images and video sequences, and marking two images. The generic head 
model without texture is used as a guide. Recorded instructions are lip-synced 
with the head directing the user to first look at a dot on the screen and push a key 
to take a picture. A second dot appears and the user is asked to take the second 
still image. The synthetic face mimics the actions the user is to follow. After the 
two still images are taken, the guide directs the user to slowly turn his/her head to 
record the video sequences. Finally, the guide places red dots on her own face and 
directs the user to do the same on the two still images. The collected images and 
markings are then processed and a minute or two later they have a synthetic head 
that resembles them. 
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Animation 

Having obtained the 3D textured face model, the user can immediately 
animate the model with the application of facial expressions including frowns, 
smiles, mouth open, etc. 

To accomplish this we have defined a set of vectors, which we call 
posemes. Like the metric vectors described previously, posemes are a collection of 
artist-designed displacements. We can apply these displacements to any face as 
long as it has the same topology as the neutral face. Posemes are collected in a 
library of actions and expressions. 

The idle motions of the head and eyeballs are generated using Perlin's noise 
functions [24, 25]. 

Conclusions 

We have developed a system to construct textured 3D face models from 
video sequences with minimal user intervention. A new head motion estimation 
algorithm takes advantage of the physical properties of human face features. The 
algorithm significantly improves the robustness over traditional motion estimation 
methodologies. It can be applied to human face modeling and tracking systems 
where the markers can be obtained either through user intervention or by using 
automatic feature detection algorithms. This algorithm can be easily extended to 
general cases where the number of feature points is not necessarily five. 

Although details of specific implementations and embodiments are 
described above, such details are intended to satisfy statutory disclosure 
obligations rather than to limit the scope of the following claims. Thus, the 
invention as defined by the claims is not limited to the specific features described 
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above. Rather, the invention is claimed in any of its forms or modifications that 
fall within the proper scope of the appended claims, appropriately interpreted in 
accordance with the doctrine of equivalents. 
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